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Jh Abstract 
m 

Recent results from Higgs boson and supersymmetry searches at the Large Hadron Collider pro- 
vide strong new motivations for supersymmetric theories with heavy superpartners. We reconsider 
focus point supersymmetry (FP SUSY), in which all squarks and sleptons may have multi-TeV 
masses without introducing fine-tuning in the weak scale with respect to variations in the funda- 
mental SUSY-breaking parameters. We examine both FP SUSY and its familiar special case, the 
FP region of mSUGRA/CMSSM, and show that they are beautifully consistent with all particle, 
astroparticle, and cosmological data, including Higgs boson mass limits, null results from SUSY 
searches, electric dipole moments, b — > 57, B s — > , the thermal relic density of neutralinos, 

and dark matter searches. The observed deviation of the muon's anomalous magnetic moment 
from its standard model value may also be explained in FP SUSY, although not in the FP region 
of mSUGRA/CMSSM. In light of recent data, we advocate refined searches for FP SUSY and 
related scenarios with heavy squarks and sleptons, and we present a simplified parameter space to 
aid such analyses. 

PACS numbers: 12.60.Jv, 11.30.Pb, 95.35.+d 
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I. INTRODUCTION 



Since its discovery decades ago, supersymmetry (SUSY) has attracted more attention 
than any other principle for physics beyond the standard model (SM). Of particular interest 
is weak-scale SUSY, which holds the promise of providing natural resolutions to the gauge 
hierarchy and dark matter problems. For the last year, the Large Hadron Collier has been 
colliding protons with protons at a center-of-mass energy of 7 TeV. The ATLAS and CMS 
experiments have each analyzed over 1 fb _1 of data and collected over 5 fb" 1 , but have not 
reported evidence for new physics p~H6]. These null results have excluded generic SUSY 
models with light superpartners and large missing Et signatures. 

Although these LHC results have disappointed the most optimistic SUSY enthusiasts, 
they do not remove the possibility that weak-scale SUSY is realized in nature. Rather, they 
shift attention to supersymmetric models that have heavier superpartners or less distinc- 
tive signatures. The former possibility is particularly natural to consider, since stringent 
constraints on flavor- and CP-violation have long motivated heavy squarks and sleptons 
of the first two generations, and experimental bounds on the Higgs boson mass have long 
motivated heavy third generation squarks to raise the Higgs boson mass through large ra- 
diative corrections. This possibility has now received even greater motivation from recent 
results from the ATLAS and CMS experiments, which combined confine the possibility of 
a light Higgs boson to the mass window 115.5 GeV < rrih < 127 GeV, and indicate excess 
events consistent with the production of Higgs bosons with masses of 126 GeV and 124 
GeV, respectively [TJ [8] . Of course, the possibility of multi-TeV third generation squarks is 
generically in tension with the requirement that SUSY resolve the gauge hierarchy problem. 

In this study, we consider focus point (FP) SUSY [9HT2] in light of recent results. We 
are motivated to consider FP SUSY for several reasons. First, in FP SUSY, all squarks and 
sleptons may be multi-TeV without increasing the fine-tuning in the weak scale with respect 
to variations in the fundamental SUSY-breaking parameters. Naturalness is a notoriously 
brittle and subjective criterion, but in this sense, FP SUSY is the unique framework that 
naturally accommodates multi-TeV top and bottom squarks. Second, many observables, 
including those at colliders, in low-energy probes, and those related to dark matter, are 
insensitive to the details of the heavy scalar spectrum, since the scalars decouple. For these 
observables, FP SUSY may be viewed as an effective theory that captures the essential fea- 
tures of a large class of models with heavy superpartners. And last, a special case of FP 
SUSY is realized in the FP region of minimal supergravity (mSUGRA) or the constrained 
minimal supersymmetric standard model (CMSSM), heretofore referred to as the "FP re- 
gion." Given the amount of work devoted to this model, FP SUSY is a practical and natural 
starting place for considering SUSY models with heavy superpartners that are newly mo- 
tivated by LHC data. For other recent work on FP SUSY and the related framework of 
hyperbolic branch SUSY [13] motivated by recent results, see Refs. [T4"t [T5]. 

We begin in Sec. |TT] by reviewing the general framework of FP SUSY and its well- 



known special case, the FP region. In Sec. Ill, we show Higgs mass predictions in 



mSUGRA/CMSSM, determine the parameter space favored by Higgs mass bounds, and 



find that current limits favor the FP region. In Sec. |IV| we show that constraints on the 
electron and neutron electric dipole moments (EDMs) are naturally satisfied in FP SUSY. 
In Sec. |Vj we then focus on the part of the FP region that has the correct neutralino thermal 
relic density Q x . This is typically presented as a thin strip in the (mo, M ly / 2 ) plane with fixed 
tan/3. To allow a more comprehensive presentation of FP results, we instead fix mo to give 
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the correct Q 



A ' 



0.23. In Sec. VI 



and in Sec. VII 



and present results in the (tan/3, M 1( / 2 ) plane, with every point satisfying 
we present results for b — > 37 and B s — » ii + ir in the (ta.11 f3, Mi/ 2 ) 
we analyze implications for dark matter direct detection and show 



plane, 

that FP SUSY remains consistent with current null results. Finally, in Sec. |VIII we show 
that the observed deviations of (g — 2)^ from SM expectations may be easily explained in FP 
SUSY (but not in the FP region). Our findings are summarized in Sec. IX The robustness 
of our numerical analyses is discussed in the Appendix. 



II. FOCUS POINT SUPERSYMMETRY 



In SUSY, the Z boson mass is determined at tree-level by the relation 
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where fi is the Higgsino mass parameter, 



m 



Ha 



are the soft SUSY-breaking Higgs mass 



parameters, tan/3 = (H®) / (H%) is the ratio of Higgs boson vacuum expectation values, and 
all of these are evaluated at a renomalization group scale near m wea k ~ 100 GeV — 1 TeV. 
For the moderate and large values of tan /3 required by current Higgs mass bounds, this may 
be simplified to 
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The weak-scale parameter m 2 H depends on a set of fundamental parameters {ai}, typically 
taken to be grand unifed theory (GUT)-scale soft SUSY-breaking parameters, such as scalar 
masses mj, gaugino masses Mi, and trilinear scalar couplings Ai. Naturalness requires that 
mz not be unusually sensitive to variations in the fundamental parameters dj. This does 
not necessarily imply ~ mz for every i, however, because terms involving some a% in the 
expression for m? z may be suppressed by small numerical coefficients. 

In the class of FP SUSY models studied in Refs. [HHI2], the fundamental GUT-scale 
parameters satisfy 
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all other scalar masses < O(10 TeV) 



Mi, Ai < 1 TeV 
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(4) 
(5) 
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for high tan /3, where the top and bottom Yukawa couplings are comparable. In Eqs. (J3]) 
and ([6]), x and x' are arbitrary constants, but for any values of x and x', the weak-scale is 
insensitive to variations in mo, even for multi-TeV mo- In other words, with these GUT- 
scale boundary conditions, renormalization group evolution takes m\ to values around m 2 z 
at the weak scale, almost independent of its initial GUT-scale value. This "focusing" of 
renormalization group trajectories does not apply to the top and bottom squark masses or, 
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of course, to any other squark and slepton masses. As a result, in FP SUSY, all squarks 
and sleptons may have multi-TeV masses without introducing fine-tuning in the electroweak 
scale with respect to variations in the fundamental soft SUSY-break parameters. For an 
extended discussion of naturalness in FP SUSY, see Ref. [TT] . 

As evident from Eqs. ^ and ([6]), the framework of FP SUSY is quite general. If one 
assumes that x — x' — 0, that all other sfermion masses are also unified to the same 
m , that all gaugino masses are unified, and that all A-parameters are unified, FP SUSY 
parameter space intersects the mSUGRA/CMSSM parameter space in what is known as the 
FP region. In general, however, FP SUSY requires neither gaugino mass nor A-parameter 
unification, and also does not constrain scalar masses that are only weakly coupled to the 
Higgs sector, such as the first and second generation squark and slepton masses. In much of 
the analysis below, we will consider the FP region, as in many cases, it serves as an adequate 
representative of general FP SUSY. The distinction between FP SUSY and the FP region 



will be relevant, however, when we discuss FP SUSY predictions for (g — 2) M in Sec. VIII 



III. HIGGS BOSON MASS 

As is well-known, current bounds from LEP2 require the Higgs boson mass to be m^ > 

114.4 GeV |16j . In SUSY, where the limit m^ < rriz applies at tree-level, large radiative cor- 
rections from heavy top and bottom squarks are required to satisfy this bound. A significant 
phenomenological advantage of the FP SUSY framework is that it naturally accommodates 
heavy third generation squarks, and with them, relatively heavy Higgs bosons consistent 
with the LEP2 bound. Given recent Higgs boson results from the LHC [7J [8], it is, of 
course, also interesting to investigate whether Higgs boson masses in the allowed window 

115.5 GeV < rrih < 127 GeV are possible, and whether masses as large as ~ 125 GeV may 
be naturally accommodated. 

In Fig. [T| we plot contours of constant Higgs boson mass in the (mo, M1/2) plane 
of mSUGRA/CMSSM. Also shown is the contour on which the neutralino relic density 
satisfies Q x ~ 0.23. Here and throughout we use SOFTSUSY 3.1.7 [T7j to generate the 
SUSY spectrum, and MicrOMEGAs 2.4 [18] to calculate the relic density and several other 
observables. In each case, we use a top quark mass of m t = 173.1 GeV and strong coupling 
constant a s (M z ) = 0.1172. 

Restricting attention to the cosmologically favored contour with Q x ~ 0.23, we see that 
the Higgs mass bound rrih > 114.4 GeV requires either mo > 2 TeV (the FP region), or 
very low mo and M1/2 ^ 500 GeV (the co-annihilation region). For the parameters plotted, 
then, the LEP2 Higgs mass bound has already eliminated much of the parameter space now 
excluded by null results from LHC SUSY searches. In the FP region, the Higgs boson mass 
satisfies m h > 114 GeV, and extends up to 122 GeV (124 GeV) for M 1/2 ~ 1 TeV (2 TeV). 
Given an estimated 2-3 GeV uncertainty in the Higgs boson mass calculation [T71 HU 120] , the 
FP region beautifully predicts Higgs boson masses in the currently allowed range from 115.5 
GeV to 127 GeV, and also naturally accommodates the 124-126 GeV mass range tentatively 
indicated by LHC search results. Varying A within the range \A \ < TeV can also raise the 
Higgs boson mass slightly by ~ 1 GeV. 

Contours of constant dark matter mass m x are also shown. Note that m x ~ (9(100 GeV), 
even for multi-TeV mo in the cosmologically-favored regions. The viable FP region contains 
heavy sleptons and squarks, but potentially sub- TeV gluinos, electroweak gauginos and 
Higgsinos as light as 200 GeV, and neutralino dark matter as light as 100 GeV, even under 
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(a) rrih and m x in GeV for tan/3 = 10 
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FIG. 1. Contours of the light Higgs boson mass m/j in black (dotted) and lightest neutralino mass 
m x in red (dashed) in the (mo, Mi/ 2 ) plane for tan (3 = 10 (left) and 50 (right), = 0, and /u > 0. 
On the blue (solid) lines, the neutralino relic density is Vt x ~ 0.23. 



the restrictive assumption of gaugino mass unification. We will return to the cosmological 



implications of FP SUSY in Sec. VII 



IV. ELECTRIC DIPOLE MOMENTS 



FP SUSY is also motivated by constraints from EDMs. Generic SUSY theories with weak- 
scale superpartners violate low-energy flavor- and CP-violation constraints. Although there 
are well-known mechanisms to suppress flavor violation, these do not typically suppress 
CP violation. In general, all gaugino masses, A-terms, and the n parameter can possess 
phases that give rise to CP violation. The most limiting CP- violating, but flavor-conserving, 
observables are the EDMs of the electron and neutron, which can arise from loop diagrams 
with either left-right sfermion mixing or a gaugino-Higgsino flip within the loop. Even with 
A 0, left-right mixing for first generation sfermions is typically negligible, but an EDM 
contribution can still arise if there is a mismatch between the phases of the gaugino masses 
and the phase of //. 

To examine these effects, we consider a simple extension of mSUGRA/CMSSM where the 
gaugino masses and fi have general CP-violating phases and the mismatch is parameterized 
as 0cp- The dominant diagrams involve left-handed sfermions and charginos with a Wino- 
Higgsino mixture, leading to contributions [21] 



df = -errif g\ |M 2 /i| tan (3 sin0cp Kc{m 2 ; , |/^| 2 , |M 2 



(9) 



where Kc is a kinematic function [22] . Diagrams involving sfermions and neutralinos produce 
sub-dominant contributions. 
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FIG. 2. Upper limits on sin^cp from neutron EDM constraints in black (dotted) and electron 
EDM constraints in green (dashed) for tan/3 = 10 (left) and 50 (right), ^4 = 0, and [i > 0. On the 
blue (solid) line, the neutralino relic density is fi x ~ 0.23. 



The current bounds on the electron and neutron EDMs are d e < 1.6 x 1CU 27 e cm [23] 
and d n < 2.9 x 1CU 26 e cm [23]. Assuming m u = 3 MeV, = 5 MeV, the naive quark 
model relation d n = (4d<i — d u )/3, and neglecting cancellations between different diagrams, 
we may derive bounds on the phase mismatch 4>cp- 

Figure [2] shows the upper limits on sin0cp in the (mo, M1/2) plane from electron and 
neutron EDM constraints. In mSUGRA, m^ r < m,7 r ~ m d L -> anc ^ so the electron EDM 



provides the stronger bound, but the neutron EDM bound is also stringent. From Fig. 2 (a 
for example, we see that for tan/3 = 10, the constraints Q x ~ 0.23 and sin^cp ^ 0.01 can 
only be satisfied in the FP region, and at the same time, the FP region with Mi/ 2 ^ 1 TeV 
can accommodate natural values of sin0cp ~ 0.3. The EDM bounds become even stronger 
for large tan/3, but may be satisfied in the FP region for M1/2 ~ 2 TeV for sin0cp ~ 0.1. 
Absent a compelling mechanism for suppressing flavor-conserving CP violation, bounds from 
electron and neutron EDMs have long ago restricted mSUGRA/CMSSM parameter space 
to the FP region, irrespective of recent LHC results from SUSY and Higgs boson searches. 



V. FP SUSY WITH FIXED RELIC DENSITY 

Results for the mSUGRA/CMSSM framework are conventionally presented as in Figs. [T] 
and [2j In these figures, the cosmologically desirable region with Q x ~ 0.23 is just a thin 
strip running through the plane, and the cosmologically desirable FP region is just a small 
part of that. Given that much of the rest of the cosmologically favored mSUGRA parameter 
space is now excluded, however, as well as our focus on FP SUSY in this study, it is more 
appropriate to consider a parameter space in which every point is in the cosmologically 
favored part of the FP region. 
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FIG. 3. Contours of (a) mo (in TeV) and (b) \x (in GeV) in the (tan/3, .M1/2) plane. Every point 
in the parameter space is in the FP region and satisfies £l x ~ 0.23, A§ = 0, and fx > 0. 



For a neutralino LSP in the FP region, a significant Bino-Higgsino mixture is required 
to produce a sufficiently low relic density, with the Higgsino component increasing with 
mo- Thus the value of mo satisfying Q x ~ 0.23 for a particular set of other parameters 
represents a lower bound. If the neutralino composes only a fraction of the relic density, 
Q x < 0.23, scalar masses are increased somewhat and the primary effect on our conclusion 
is a weakening of direct detection limits. It is also possible to disconnect the FP effect on 
fine-tuning from cosmological considerations by introducing a gravitino LSP which allows 
a larger neutralino relic density to be considered; we restrict our intention to the case of a 
neutralino LSP. 

To satisfy the relic density constraint, we continue to consider fixed values of A Q and 
sign(/i), but require the neutralino to be a thermal relic with Q x = 0.23. This implies a 
constraint on the remaining parameters mo, M1/2, and tan/3. We choose M1/2 and tan/3 as 
the free parameters, and use VL X to determine m-o[j] In general there are several values of mo 
satisfying this condition for a particular (M1/2, tan /?) pair, arising from the co-annihilation 
region at low mo, the FP region at large mo, and the A- funnel region for moderate mo and 
large tan 0. We focus on the FP region by always choosing the largest value of m for a 
given point in the (Mi/ 2 ,tan/3) plane. 

In Fig. [3j we show contours of constant m and fx in the (tan (3, M 1 / 2 ) parameter space 
defined above, where every point has Q x 



we see 



0.23, A = and fx > 0. In Fig. [3(a) 
that mo increases as Mi/ 2 increases and decrease as tan/3 increases. In the FP region the 
large mass of the sfermions makes them nearly decoupled for the relic density calculation. 
The correct value of mo is instead solely determined by its impact on the Higgs potential, 
which sets \ft\, and which in turn determines the correct Higgsino-Bino mixture to produce 
0.23n In Fig. 3(b) , we see that fx grows with increasing M1/2, but is nearly independent 



1 Alternatively, one could choose Mi/ 2 a nd mo as the input parameters, and predict tan/3 |25j . 

2 The determination of too in Fig. 3(a) is sensitive to the value of the top mass (see, e.g., Ref. [26]), and 
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tan (3 tanp 

(a) m h (GeV) (b) m x (GeV) 

FIG. 4. Contours of (left) and m x (right) in the (tan/3, M1/2) plane. Every point in the 
parameter space is in the FP region and satisfies £l x ~ 0.23, Aq = 0, and n > 0. 

of tan (3, given the subdominance of terms involving tan (3 in the neutralino mass matrix. 

In Fig. |1J we plot contours of mn and m x in the same (tan (3, M 1 / 2 ) parameter space. 
The large value of in the FP region raises the Higgs mass well above the LEP2 bound of 
114.4 GeV, and is confined to the currently allowed range of 115.5 GeV < < 127 GeV. As 
one moves to smaller values of tan 0, mh increases even though its tree-level value drops, be- 
cause of the enhancement of the loop-level contribution from increasing m . The neutralino 
mass contours satisfy m x ~ M\ ~ 0.4Mi/2, since the neutralino is primarily Bino-like, al- 
though there is an increasingly significant Higgsino component as M1/2 increases. As with 
the Mi and fi contours, the m x contours are also nearly independent of tan /3. 



VI. RARE B PROCESSES 

Rare decays are often used to constrain new physics scenarios, and in particular, the 
decays B — > Xgj and B s — > are well-known probes of new physics. The measured 

value of B(B A s7 ) is (3.55 ± 0.33) x 10" 4 [28], consistent with the SM value of (3.15 ± 
0.23) x 10~ 4 [SU [30]. The value of B(B S — > /i + /i _ ) has been the subject of recent interest, 
with a CDF analysis reporting a value of 1.8^o'g x 10~ 8 , and claiming 4.6 x 10~ 9 < B(B S — > 
fi + fi-) < 3.9 x 10~ 8 at 90% C.L. [311 [32]- Meanwhile, CMS and LHCb analyses produced 
only upper limits at 90% C.L. of 1.9 x 10 -8 [33] and 5.6 x 10 -8 [M], respectively, and 
1.08 x 10~ 8 [33] from a combined analysis using 2010 LHCb data [36J. The SM value is 
(3.19 ± 0.19) x 10~ 9 [33 EH], consistent with the LHC bounds and marginally inconsistent 
with the CDF analysis. 

varies somewhat for different MSSM spectrum generation programs. The determination of /i shown in 
Fig. |3(b)[ however, is preformed directly from a fit to the measured relic density and is thus robust and 
independent of the value for the top mass or the spectrum generator used [27]. For more details, see the 
Appendix. 
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FIG. 5. Contours of AB(b -> s 7 ) in units of 10~ 4 (left) and AB(B S -> ^V") in units of 1(T 8 
(right) due to SUSY in the (tan/3, M^) plane. Every point in the parameter space is in the FP 
region and satisfies O x ~ 0.23, Aq = 0, and \i > 0. 



Figure [5] shows the contributions to B — > X s ^ and B s — > [i + [iT from supersymmetric 
particles. For both observables, the primary supersymmetric contributions arise from loop 
diagrams involving either charginos or charged Higgs bosons. For B — > X s ~f the former 
produces a suppression in the decay for ji > and an enhancement for /i < and the latter 
an enhancement for either sign of fi. For B s — > fi + fi~ the chargino contribution is negative 
and charged Higgs contribution positive for either sign of /i. Within the FP region the 
chargino diagram dominates. For B — > X s 7, this puts the supersymmetric result in greater 
tension with experiment than the SM result for fi > 0, though only significantly so at low 



M1/2 and large tan/3 — the 2a discrepancy line is plotted in Figure 5(a) For fi < the 



contribution is positive and within la of the observed result for the entire [M1/2, tan /3) plane. 
For B s /i + yU _ , the supersymmetric contribution in the FP region does not significantly 
alter the SM prediction, at least relative to current experimental uncertainties. 



VII. DIRECT DETECTION OF DARK MATTER 



In the cosmologically-favored region of the FP region, neutralinos make up the dark 
matter. These regions of parameter space are then constrained by null results from dark 
matter searches. In particular, null results from direct detection searches that constrain the 
spin- independent x- nuc l eon cross section a p have been advanced as significant constraints 
on FP SUSY [MSI]. 

In the absence of large left-right mixing, the dominant contributions to both neutralino 
annihilation and spin-independent scattering are dependent on the "Higgsino-ness" of the 



9 



> 

CD 




ix0.5 







-°* 


-in 

- 2 DM 





-In 

5 DM 



100 200 300 400 500 600 700 800 900 100C 



rr\ (GeV) 



(a) Higgsino content ( 10 ) for various relic densities 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
fs 

(b) Dependence of a p on the strange quark form 
factor f s 



FIG. 6. Left: Higgsino- ness ag for neutralinos in the FP region consistent with various relic 
densities. The right-most point on each curve corresponds to M 1 / 2 = 1 TeV. Right: The spin- 
independent X" nuc l eon cross section function of f s for a model in the FP region with 
(too, Mi/ 2) = (3 TeV, 550 GeV). The shaded regions indicate the la uncertainties on the various 
f s determinations. In both plots, tan /3 = 10, Aq = 0, and fx > 0. 



lightest neutralino, defined as 



"H 



at, 



+ a 



where the neutralino eigenstate is 



X = a s B + a w W + afj u H u + a^Rd 



(10) 



(11) 



Figure 6(a) shows the dependence of ag on m x in the FP region. The Higgsino-ness 
generically increases with m x to offset the suppression in annihilation from the lowered 
cross section. However, it decreases when new annihilation channels open at m x ~ mw,fnz 



and m x ~ m t . Figure 6(a) also shows curves in which the neutralino makes up only a fraction 



for lower relic densities, ag increases to enhance the annihilation rate. 



10, A = 0, and 



of the relic density 

The curves are generated by varying My 2 up to 1 TeV, for fixed tan/3 
/j > 0. 

To determine the spin-independent x-nucleon cross section a p , the contributions of the 
couplings to each individual quark must be considered. The individual couplings must be 
weighted according to the scalar quark form factors f^, typically parameterized as 



[N\m q ip q ip q \N) 



f a N M N 



(12) 



The parameters f^ d are reasonably well known, and the heavy quark contributions are set 
by loop contributions using the gluon form factor. However, the value of f^ is less certain, 
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given discrepancies between current experimental and lattice results, and this is a well- 
known source of uncertainty for direct detection predictions [4*0] W2\ "4"3] . The experimental 
determination combines a derivation of the pion-nucleon sigma term from meson scattering 
data jSj combined with a number of chiral perturbation theory results [loT - PTT] . giving 

/. = /; = fS ~ 0.36 . (13) 



More recent calculations support older determinations of the pion-nucleon sigma term 
For this value of f s , the other form factors are all much smaller, f£L s < 0.05, and so the 
strange quark contribution dominates the direct detection cross section [12]. However, two 
recent lattice studies have found much smaller values for f s [49J 150], with an average of 
f s ~ 0.05. For this value of f s , the strange quark contribution is much closer to that of the 
other quark flavors [43l [51] . 



Figure 6(b) shows the dependence of a p on f s for both positive and negative /x in the 
FP region. The value of a p varies by a factor of ~ 3 between the experimental and lattice 
determinations of f s , which has significant implications for direct detections bounds. The 
scattering cross section may also be suppressed if ji < 0. This possibility is often ignored 
in studies that assume /x > to reduce the discrepancy in (g — 2) M between the SM and 
experimental data. 

Figure [7] shows contours of a p for positive and negative /x and f s consistent with exper- 
imental and lattice results. The general factor of ~ 3 due to different values of f s is once 
again apparent. The cross section a p for /x < also shows a general suppression relative to 
that for /x > 0, though the suppression varies significantly with both mass scale and tan/?. 
For /x > there is a general enhancement in a p at low tan (3 due to the coupling to the light 
Higgs, and at high tan (3 due to a reduction in the masses of the heavy Higgs bosons. These 
effects are also present for /x < 0, but instead produce a suppression. 

Figure \f\ also shows the regions of parameter space excluded by XENON100 [52]. For 



/x > 0, fs = .05 in Fig.|7(a) M 1/2 < 450 GeV is excluded for all tan/3. The case of f s = 0.36 
in Fig. 7(b) is markedly different, with exclusion up to M1/2 ~ 1.6 TeV for all tan (3 and 



larger M 1 / 2 for low and high tan (3. The same trend carries over to /x < — in Fig. 7(c) the 
exclusion is limited to a small range of M1/2 at moderate tan/3 where peaks, and to a 
small region at high tan where scattering is dominated by the heavy-Higgs boson mediated 



process. The exclusion in Fig. 7(d) is greater due to larger f s but still reduced compared to 
the /x > case. 

In summary, we find that FP SUSY is far from excluded by current direct detection 
bounds. For large f s ~ 0.36 and /x > 0, significant portions of the parameter space are 
excluded, requiring M1/2 ^ 1.6 TeV to survive. For the smaller values of f s favored by 
lattice results or /x < 0, a much larger portion of the parameter space is viable, including 
regions with gaugino masses as low as M1/2 ~ 250 GeV. At the same time, it is, of course, 
interesting that the direct detection bounds are within factors of a few from probing all of FP 
SUSY. To the extent that LHC SUSY and Higgs boson results motivate SUSY with heavy 
squarks and sleptons, they also motivate direct detection experiments that are approaching 
sensitivities to zeptobarn spin-independent cross sections in the near future. 



VIII. THE ANOMALOUS MAGNETIC MOMENT OF THE MUON 

The well-known ~ 3a discrepancy between the experimental and SM values in the anoma- 
lous magnetic moment of the muon [53H55] is currently among the most compelling pieces of 
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FIG. 7. Contours of o v in zeptobarns for /i > (top) and \x < (bottom), with / s = 0.05 (left) 
and f s = 0.36 (right). In each panel, the shaded region is excluded by XENON100 [52] . 



evidence for new physics. The supersymmetric contribution is given by jl — x° and — x 
loop diagrams. The (g — 2) M discrepancy has two robust implications for SUSY — it is the 
primary result motivating relatively light superpartners, and it favors fi > 0. 

The large sfermion masses in the FP region produce too small a value for A (g — 2)^ USY 
to explain the observed discrepancy of (2.9 ± 0.9) x 10~ 9 |54j . Figure [8] shows the value of 
A(g — 2)^ USY in the FP region parameter space. The largest value attained is A(g — 2)^ USY « 
0.5 x 10 -9 , insufficient to produce even 2a agreement with the experimental result. 

As noted in Sec. [TTJ, however, FP SUSY is far more general than the FP region. In 
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tanp 

FIG. 8. Contours of the supersymmetric contribution to (g — 2)^ in units of 10 -9 . Every point in 
the parameter space is in the FP region and satisfies J7 X ~ 0.23, Aq = 0, and \i > 0. 

particular, in FP SUSY, the smuon and muon sneutrino need not have masses unified with 
the other scalars, and so may be much lighter than the third generation squarks. To explore 
this possibility and its implications for (g — 2) M , we consider the slight modification of 
mSUGRA/CMSSM in which all scalars have GUT-scale mass mo, except for the smuons 
and muon sneutrino. This modification is intended to be schematic, demonstrating the 
behavior of (g — 2) M with lowered smuon masses without bias toward a particular approach. 
A fully consistent approach must consider flavor and GUT unification issues. For simplicity, 
we take the smuon masses to be degenerate at the weak-scale, with physical masses 

M A = m H = m~ m = m p „ . (14) 

At each point in the (M l7 / 2 , tan/3) plane, we determine the value of Mp, that gives A(g — 
2)^ USY that either brings the theoretical prediction into complete agreement with the central 
experimental value or reduces the discrepancy to 2a. Note that the dominant factor in the 
determination of the relic density is the Higgs soft mass, with the sfermion masses providing 
subleading effects, as long as rriq > 500 GeV and > 200 GeV [56J. The smuons can 
therefore be quite light without affecting the relic density constraint. 

The results are given in Fig. |9j As M\/2 increases, the required smuon mass decreases 
to maintain a constant SUSY contribution to (g — 2) M , and at some point, the required Mp 
becomes too low, as it implies a /2 LSPrl The supersymmetric contribution A(g — 2)^ USY 
also has a linear dependence on tan (3, and so at large tan (3, there are allowed solutions for 
larger values of M 1 / 2 and M^. 

It is important to check that the scenarios for resolving the (g — 2)^ discrepancy are viable 
in light of null results from LHC new physics searches. The model-independent bounds on 
slepton masses are, of course, far weaker than those on squark masses. The best limits on 
slepton masses are still those from LEP2, which require > 100 GeV [57]. In the future, 



3 Mp, is cut off at 1.1 x ra x numerically to avoid recalculating the relic density due to /2 — x° co-annihilation. 
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FIG. 9. Contours of required to eliminate the discrepancy between the theoretical and experi- 
mental values for (g — 2)^ (left) and to reduce the discrepancy to 2a (right). This model framework 
is a slight modification of mSUGRA/CMSSM in which all scalars have GUT-scale mass mo, except 
for the smuons and muon sneutrino, which have physical mass Mn. In the shaded regions, the 
fl becomes the LSP. To specify all parameters aside from the smuon and muon sneutrino masses, 
every point in the parameter space is in the FP region and satisfies fi x ~ 0.23, Aq = 0, and \x > 0. 



with 30 fbT 1 of data at 14 TeV, the LHC will be able to discover sleptons through Drell- 
Yan production for m^ L < 300 GeV and < 200 GeV [SB]- Greater sensitivity may be 
available in scenarios where the sleptons are produced in cascades jSH] ■ However, in the FP 
region where all other scalars are heavy and gluino production dominates, if the sleptons are 
heavier than all charginos and neutralinos, they will not be produced in gluino cascades, and 
so the Drell-Yan limits apply. This is the case for regions of the (tan /3, Mi/ 2 ) plane shown in 
Fig. |9j and so there are viable FP SUSY scenarios that resolve the (g — 2) M discrepancy. It 



would, however, be interesting to investigate scenarios motivated by the (g — 2)^ discrepancy 
in which sleptons are produced in gluino cascades. 



IX. CONCLUSIONS 

SUSY models with heavy squarks and sleptons have long been motivated by constraints 
on flavor- and CP-violation, the LEP2 constraint on the Higgs boson mass, and other con- 
straints, such as proton decay bounds. Recent null results from LHC SUSY searches have 
further focused attention on this possibility, and the interest in such scenarios is especially 
heightened by the currently allowed Higgs boson mass window 115.5 GeV < < 127 GeV, 
and tentative indications from the ATLAS and CMS experiments for a Higgs boson with 
mass near 125 GeV. 

Generic SUSY scenarios with heavy sfermions, and particularly heavy top and bottom 
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squarks, imply fine-tuning of the weak scale, subverting the basic motivation for weak-scale 
SUSY. In FP SUSY, however, this is not the case. The mass parameter m 2 H evolves to 
values around m 2 z at the weak scale, almost independent of its GUT-scale starting value. 
This focusing of RG trajectories implies that the weak scale in FP SUSY theories is not 
fine-tuned with respect to variations in the fundamental SUSY-breaking parameters. Note 
that the fact that m 2 H evolves to values around m 2 z at the weak scale for a particular choice 
of GUT-scale parameters is necessary to remove fine-tuning with respect to variations in /i, 
and is possible for other choices of GUT-scale parameters (see, for example, Ref. [601 161]). 
However, naturalness with respect to variations in all SUSY-breaking parameters requires 
that rn 2 Hu evolve to a weak-scale value irrespective of its starting value, and so the focus 
point behavior of renormalization group trajectories is an essential feature of any natural 
theory with multi-TeV top and bottom squarks motivated by the currently allowed Higgs 
boson mass range. 

In this study, we have focused for the most part on models of FP SUSY that are also 
part of the mSUGRA/CMSSM framework. These FP region models naturally produce Higgs 
boson masses above the LEP2 bound of 114.4 GeV, and suppress electron and neutron EDMs 
sufficiently, even for 0(1) phases. To more globally display the predictions of FP SUSY, 
we have required fl x ~ 0.23 and plotted results in the (tan/3,Mi/ 2 ) plane. We find that 
FP SUSY naturally accommodates Higgs boson masses up to 120-124 GeV, which, given 
an estimated 2 GeV uncertainty in the theoretical calculation, is consistent with current 
Higgs boson mass indications. In addition, we have shown that FP SUSY is naturally 
consistent with constraints from b — > S'j, B s — > n + fi~, and null results from dark matter 
direct detection experiments. Finally, in general FP SUSY with a non-unified smuon mass, 
we have found that FP SUSY may resolve the discrepancy in (g — 2)^ consistent with all 
current constraints. 

Given these successes, it is natural to ask what evidence for FP SUSY should accumulate 
in the near future if FP SUSY is realized in nature. Certainly the Higgs boson should be 
discovered with a mass in the currently allowed mass window, and searches for SUSY from 
gluino pair production, followed by gluinos cascading through charginos and neutralinos are 
promising for some of the parameter space [62H66] . Equally exciting would be the discovery 
of dark matter with a spin-independent x-nucleon cross section near the zeptobarn scale, 
which is a robust prediction of mixed Higgsino-Bino dark matter with heavy squarks and 
sleptons. Finally, most signals of indirect dark matter detection are also generically enhanced 
in the FP SUSY scenario [67J. 
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APPENDIX 



It is well known that different spectrum calculators do not give identical results for the 
SUSY mass spectrum, even for the same set of input parameters [68l [69] . The reasons for 
this apparent discrepancy are well understood; see Ref. [27] for a nice summary. Above all, 
one should keep in mind that the SUSY spectrum is always calculated at a fixed order in 
perturbation theory, and there is an intrinsic uncertainty due to neglecting the higher order 
terms in perturbation theory. The main differences between the various programs arise 
mostly because they choose to neglect different sets of higher-order terms. For example, one 
may choose to use either tree-level or 1-loop-corrected masses in the radiative corrections, 
or choose a slightly different value for the matching scale between the SM and the MSSM. 
In each case, the difference between the two options is a higher-order effect. In this paper, 
we chose to work with the SOFTSUSY program, but we expect qualitatively similar results 
from other spectrum generators as well. 

On a related topic, each spectrum calculator needs to solve a two-sided boundary value 
problem, since the boundary conditions for the gauge and Yukawa couplings are specified 
at the weak scale, while the soft SUSY-breaking parameters are given at the (yet to be 
determined) GUT scale. The standard approach used by all programs is to apply iterations 
until converging on a solution. Unfortunately, on occasion one may encounter poor conver- 
gence as a sign of a chaotic behavior 
through the (m ,Mi/ 2 ) plane of Fig. 2(a 



Oj. This is illustrated in Fig. 10, which takes a slice 



850 GeV. The figure shows the chargino mass M-+ 



Of Mi/a 

(right axis) and the relic abundance calculated by MicrOMEGAs (left axis) 



in 5 GeV increments along mo, for a fixed value 

calculated by SOFTSUSY 
We see that 
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at low mo, SOFTSUSY is able to converge, and both quantities follow a well-defined trend. 
However, at sufficiently large values of m , SOFTSUSY is not able to achieve the desired 
level of convergence, and the obtained results (upon exiting after a fixed number of itera- 
tions) visibly deviate from the expected trend. As seen in Fig. 10 , in principle this presents 
a problem for the correct mapping of the boundary of the region allowed by LEP chargino 
searches (M-+ > 103 GeV). Fortunately, however, the parameter space points with the 
desired value of the relic density (Q x h 2 m 0.1) are relatively safe, since they are still well 
within the region with good convergence, and the maps shown in Figs. [3] and [4] are robust. 
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